A(nother) Continuum Model for Dephasing in Mesoscopic Systems 
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A dephasing model in the spirit of Biittiker's fictitious probe model where infinite probes are 
distributed uniformly over the conductor is proposed. The dephasing rate enters into the model as an 
adjustable parameter and to compute the conductance. A one-dimensional delta function scatterer 
model is solved numerically. We observe the dephasing effects on the calculated conductance. 

PACS numbers: 



I. INTRODUCTION 



Dephasing, the loss of coherence in wavefunction, is an important phenomenon in mesoscopic systems. It is the 
phenomenon that distinguishes the microscopic where full quantum coherence is the rule and the macroscopic where 
there is no trace left of the quantum phase. In the intermediate mesoscopic regime, its effect is important. It is 
either due to the collisions with the other electrons and phonons, which can be adjusted by temperature or it can be 
I——, influenced by external factors^. 

]~~| Several models have been proposed for modelling the dephasing effects, coherent absorption, wave attenuatioiP, 

introducing random phase fluctuations in the scattering matrix^ are a few. One of the oldest models is the flctitious 
1^ probe model of Biittiker.* ^ This model has been applied into several different problems. Also, it has been changed as 
c/3 a model to overcome some of its deflciencies; for example momentum randomization is eliminated and pure coherence 
effects are brought to froniP and long stub model is applied for satisfying the charge conservation requirement for time 
\^ depen dent c urrents.*^ (But the stub model is introduced earlier). The model can be justifled based on microscopic 
^ theory.lS12112l They are also generalized to the continuous case where infinite probes are distributed continuously over 
g the conductoiMimm. 

In this contribution, we are going to propose another model based on Biittiker's fictitious probe model where infinite 
' probes are distributed continuously over the conductor. The inelastic scatterers are modelled in terms of a scattering 
^ matrix with a coupling parameter D, which sets the strength of the decoherence introduced. The aim of this paper is 
Q to present this continuous model in order to get the conductance of a one-dimensional conductor. In the next section 
^ we define the discrete model and after that section the continuum version of it and numerical procedure is given. The 
last section is devoted to our results and conclusions. 



II. THE DISCRETE MODEL 



> 
00 

l> We are interested in extending Biittiker's model for decoherence in ID transport in a way that decoherence proceeds 

at every location. The geometry of the problem is shown in Fig. [T] Here there is a conductor along which electrons 
move and scatter. Apart from that, N additional probes are also placed for modelling the decoherence effects on 
the main conductor. It is assumed that the electrons can jump between the conductor and the probes. It can go to 
equilibrium in those probes but will eventually return back and at the end coherence with the wavefunction in the 
main conductor will be lost. 

In order to describe the possible states of the electrons, the state of electron at position x on the main conductor 
\^ is denoted by \x) and the state when the electron is on probe-j at the position ^ will be denoted by j). Any state 
I'f/') can be expressed as a superposition of these as 

1^) = j dxi;{x) \x) d^MO IC, j) , (1) 

where tpix) is the wavefunction on the main conductor and (j^ji^) is the wavefunction on probe-j. We let the potential 
on the main conductor be V{x). On the probes, we assume that the electrons move freely, feeling the constant 
potential Vj on probe-j. 

The Hamiltonian for the electrons is taken as 



H^ho + Y,h,+J2 dj (IC = 0, j) {x,\ + \x,) (e - 0, j|) 

J 3 



(2) 
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FIG. 1: The geometry of the problem. 



where ho and hj denote those Hamiltonians of the main conductor and probes respectively and dj is a real number 
representing the coupling strength to probe- j'. 

We can write down the following differential and abstract representations of ho and hj 

'^0 =-£^£2+V{x) , ho ^ J J dxdx'ho{x;x')\x){x'\ , (3) 

--&ll^ + Vj , h, = j j d^d£,'h,{t,a\^,j){^\j\ ■ (4) 

Note that these operators act on their respective spaces. As a result, we have ho I'?, j) — hj \x) = 0. Also hj i) = 
if i ^ j. As a result, for the state given in Eq. we have 

holij) ^ Jdx(^-^^^ + V{xMx)^\x) , (5) 

h' d^.iO 



^3 \^) - d^-— + (0 le, j) . (6) 



In the Hamiltonian we add a term for the transfer of electrons between probes and the main conductor. We assume 
that when the electron is at position Xj on the main conductor, it can jump to the origin, ^ = 0, of probe-j. The term 
in the Hamiltonian of the form |^ = 0,j) {xj\ handles this. The hermitian conjugate handles the opposite process, 
namely jumping from probe-j to the main conductor. 

Here dj could have been chosen complex valued, but this is unnecessary since it does not introduce any new 
effects. Moreover, the reality implies a simple time-reversal operation (complex conjugation of wavefunction) and the 
symmetry implies that the scattering matrix is symmetric. 

The Schrodinger's equation, H [ip) = E \ip) can be expressed in terms of wavefunctions as 



h^ d^i;{x) 
2m* dx"^ 

h' d^cf>,{0 



+ V{x)ij{x))+Y,djcl>j{0)S{x~x,) = E^x) , (7) 
^ j 

2 +V,ct>m)+dMxj)S{0 = Ecly.iO . (8) 



2m* dC 

where we assume that the potential on the main conductor, V{x), is constant outside a certain interval. 

V{x) ^[Vl for x< 

I Vr for X > Xj^ 

where between the points x^^^ and x^l^\ V{x) varies. The scattering region and the points Xj are contained in this 
interval. 
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Wc will define the incoming wave amplitudes aj, a'j and the outgoing wave ampUtudes bj, b'j {j = 0,1, N) for 
any solution at energy E by 

f (aoe^*^^^ + boe-'''^=') for x < xi^^ 
V'(^) = „ ,J. . , f.) (10) 



^ (a[,e-*'^«"= + ^^^e**^"^) for x> x 



R 



(n - ""^^''^ + ^^^'^"'^ for ^ < 

f'il^J - a' e-**=^« + fe^e^'^J^ for ^>0 ^ ^ 



where for any energy E, the left and right wavenumbers are defined as 



2m*{E-VL) 2m*{E-Vn) 
k, = J , kn = J . (12) 



For the probe-j, the electrons move freely with wavenumbers 



The corresponding velocities are defined accordingly, vl = fikL/m* etc. 

There are 27V + 2 independent solutions of the wave ccjuation. Any particular solution can be obtained by choosing 
arbitrary values for the incoming wave amplitudes Oj and a'^ . From these values alone, the outgoing wave amplitudes 
bj and b'^ can be determined. The relation between the outgoing and incoming amplitudes involves the scattering 
matrix 

AT 



^3 = ^S^ai + Sjvol^ , (14) 

AT 

^3 = ^^,'iai + 5,va'i . (15) 



i=0 

Our purpose is to obtain the scattering matrix. Through this we can calculate the transport properties of the system. 



A. Solution for probe-j 

First we write down the solution of the Schrodinger's equation for probe-j. The wavefunction (/>j(^) is continuous 
at the origin ^ = 0, but its derivative has a discontinuity 



A<^;.(0) = <^;.(0+)-<^;.(0-) = — ^V'(^,) . (16) 



2m* 

The outgoing amplitudes then can be expressed as 

= aj-iDj^{xj) , (17) 

bj = a'j - iDjiPixj) , (18) 

where 

Dj = -^ . (19) 

Since dj has dimensions Energy x Length, Dj has the dimensions of square root of velocity. We will need the following 
expression below. 

<^,(0) = , (20) 



where 

dj = D,i;{x,) + iiaj + a'^) . (21) 
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B. Solution for the main conductor 



Schrodinger's equation for the main conductor can be expressed as 

[E-ho]i^ix)^^dj(l)j{0)6{x~Xj)^-'ihJ2Dj0jS{x~Xj) . (22) 



which can be solved easily by using the Green function as 



i>{x) = Mx) + J dyG{x; y) i^ihY,Dje,5{y - x.Yj , (23) 

where V'o is a particular solution of the homogeneous equation, [E — /iqI'^/'o = 0; Sind G{x;y) is the Green function 
satisfying 

[E~ho{x)]G{x-y)^5{x^y) . (24) 

The wavefunction ipix) is 

^{x) = aQipL{x) + a'QipR{x) — ifi'^^G'^^\x\Xj)Dj6j . (25) 

3 

where ^l{x, E) and ^PFt,{x, E) are two scattering solutions of the main conductor. The general solution of the homo- 
geneous equation can be expressed as a superposition of these two. These solutions satisfy 

h= [e^^^'-' + roe-''^^^) for x < x^^^ 



M^,E) = ' (26) 

I 751i ^oe lOr X > Xj^' 

( t'e-'^^'-" for X < x'^j^ 

Vr{x,E) = <^ ^ ° ^ f., (27) 

These are the solutions of [/iq — -B]</?L,i? = obtained when there are no probes connected. Here rg, Tq, and t'^ 
are reflection and transmission amplitudes and we have tg — tg due to the symmetry of the scattering matrix. Green 



functions can be expressed in terms of these solutions, ^pl,r- Note that in Eq. (23 1, the term containing the Green 
function can have only outgoing waves if G'^-* is used. In that case, all incoming waves should appear in ^o- As a 
result we have V'o = ao<PL + d'o'PR- 

Since 9j depends on ipixj), we need to solve this equation. To simplify the notation we first define O^j as 

Oaj = DjM^j) + + a'j) , (28) 
and note that 9oj depends only on incoming wave amplitudes. Using this, we get the following set of N equations, 

9e = eoi~ihJ2DfG^^H^i,Xj)Djej . (29) 

Let us now define an iV x matrix T(j as 

= Stj + ifiDeG^+\xe,Xj)Dj = Sij + —!j^ipR{xj^)(pL{xj^) (30) 

to 

= Stj + -r^RK^Lj^ , (31) 
to 

where Jrj — Dj(pR{xj) and f^j = DjipL{xj). The final solution is 9j = (F^^) 6*0^ from which we obtain all 
scattering amplitudes. 

It may be shown that the inverse of T can be written as 

(r-i)^^. = % - ^ TRj^ TLj^ , (32) 

where 

TL=r-Vi , TR^T-^fR , U^to- flT-^jL . (33) 
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C. The scattering matrix 

We look at the behavior of ip{x) for x < x^^\ 



^ aoe^'"" + ^ [roao + ha'^ - > ] tmOw ) 6"*'=-^ (35) 



Therefore we have 



For X > x^j^ we get 



bo = roao + ioao - X] '^LgOoi . (36) 



6o = <oao + foflo - X tr^Om ■ (37) 



The equations ( 17|18 1 give the outgoing amphtudes at the probes as follows 

bj = -aj-i{T-')^^eoi , b'j = -a'j-i{T-')^^0oi (38) 
Finally, 6oe depends only on the incoming wave amplitudes through 

Ooe = aafu + Oo/fl^ + + ■ (39) 

From these expressions we can read off the scattering matrix elements as follows. First scattering amplitudes for 
the main conductor 

td = Slr = Srl = to - fl^^^fn. = to - tI/r ^ to - t^Jl , (40) 
Sll = r-o - /Jr- Vl - ro - , (41) 

Srr = r'a- flT-^fR ^r'^- r^fn ■ (42) 

We will use the symbol td for the amplitude Slr and call it the direct transmission amplitude. For the scattering 
into and between the probes we have 

Slj = Slj' = SjL = Sj'L -"i-TLj , (43) 

Srj — Sbj' = SjR ~ Sj'R = -iTRj , (44) 

S,,^S,>,, = -^6,,+ {T-')^^ , (45) 

V - S,,, = (r-i)^.^ . (46) 

Note that j and j' denote the negative and positive axes respectively on probe- j. These two directions are entirely 
equivalent for scattering. Therefore if an inversion is taken on probe-j (i.e., j is switch with j') then the scattering 
matrix should remain invariant. This symmetry can be seen in the expressions above. 

For example, for the scattering between two different probes j and £ (j ^ £), the scattering amplitude is 

- ^ ^Rj< ^ij> (47) 

independent of the directions the wave comes and goes. If a wave coming from probe-j is scattered back into the 
same probe (perhaps through passing to the main conductor) then the transmission amplitude is 

- 1 - r^L.TR, (48) 
td 



and the reflection amplitude is 



l + (r-'),j = -^rL,r^, (49) 
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Note also that the 5-matrix has to be unitary. An interesting question is this: Which properties should the T 
matrix satisfy so that the resultant 5'-matrix is unitary? It appears that the following equations 

(Pl{x)* = r^ifiLix) + tlipR{x) , (50) 
'PRix)* ^t*QLPL{x) + r'f^*(fBix) , (51) 

which are also satisfied by /l and fa, are the only ones we need. From here, it can be shown that F-matrix and its 
inverse satisfy 

F + rt-2/ = hfl + fnA , (52) 
r-i + (rt)-i - 2r-i(Ft)-i = Tirl + r^r^ . (53) 

Unitarity of ^'-matrix follows from these. 



D. Probabilities 



We are using mostly the transmission probabilities. The direct transmission probability is Td = \td\ ■ The trans- 
mission probability from left lead to a direction in probc-j and the corresponding quantity for the right lead are 

Tlj = \rL,? , (54) 
Tb, = |TK,f . (55) 

The transmission probabilities between two different probes j and £ can be expressed in terms of the quantities above 

2 



td 



(56) 



In other words, knowing the transmission probabilities for the main conductor, we can determine these probabilities 
between the probes. 



E. Conductance 



We suppose that the leads of the main conductor have electrostatic potentials and Wr. We assume that both 
directions on the probes are at the same potential Wj. The differences in chemical potentials are related to these 

potentials by /i^ — ^ir = {—e){WL — Wr) etc. 

The current that enters from the lead a and go to the lead (3 can be expressed as 



Ia^l3 = 2 



2e2 



^ {lic.-tJ^p) = ^{W^-W0) = Go{W^-W0) , 



(57) 



where Go is the conductance quantum. Form these we can get expressions for the total current going into a lead. 



Go 



Ir = Go 



Td{WL -Wr)+Y, 2TLjiWL - Wj) 



Td{WR -Wl)+Y^ 2Tr,{Wr - Wj) 



Ij = Ij, = Go [TLj{Wj - Wl) + TRj{Wj - Wr) + Tjj,{Wj - Wj) 



^2TejiW,-We) 



(58) 

(59) 
(60) 
(61) 



The total current going in has to be zero: II + Ir + J2j = 0. Also, all the potentials Wa can be shifted by a 
constant amount, W^ Wa + SW, and this does not change the value of currents. Due to this we can choose one of 
the potentials (such as Wr) to be (grounding). 
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Since probes are only imaginary, we require them to carry no current, Ij = 0. In this way, if electrons go into one 
of these probes, same number of electrons come back. In this way, particles do not disappear on the average on the 
main conductor. In that case we have II = —Ir = I, the current passing from the device. We suppose that Wr = 
and express all other potentials in terms of Wl- 

The equation for the current entering into probe-j is 

TljWl = Tlj + Tn, + 2 ^ T,, W,- 2 ^ T.eWe . (62) 
\ ijtj J W 

The terms inside the parentheses is equal to (by the unitarity of 5-matrix) 



1 _ \s..f - .f = ( + ( * - 2^Mk = rnj - 2^^ 



(63) 



We define a new N x N matrix, P, with 



Pje = mjdje-2'^^^i^^ . (64) 



It is a symmetric matrix with real elements which also satisfies (because of the way the diagonal elements are defined) 

^Pje = nj+Tnj . (65) 



Using this matrix, we can find the potentials W. 



Wl 

Using these, the dimensionless conductance can be expressed as 



J' 

'^'-T.Pi/Tu . (66) 



= Trf + 2 ^ - 2 ^ n,Pr/Tu (68) 
= Ta + 2^TnjP7/Tu (69) 



III. THE CONTINUUM VERSION 



We are now going to pass from the discrete model solved above to a continuum model where the probes are infinite 
in number and they are distributed uniformly to every position. Still, we may want to keep a finite range for the 
positions where these probes are in contact with the main conductor. For this reason, we suppose that the region 
where decoherence occurs is on the interval between positions and x^. 

Second, we are going to make a connection with the previous discrete problem. So, we are going to select N points 
uniformly within the decoherence interval. 

x2 < Xi < X2 < ■ ■ ■ < xn < x^ . (70) 

We are not going to specify how these points are chosen, but in A/' — > oo limit, they should fill out the whole interval. 

Let Axj be the length of interval where the point Xj corresponds to. A possible choice might be Axj = .Xj+i — Xj 
and Ax AT = x^ — xn H xi = x^- Another possibility is choosing xj in the middle of each subinterval of length Axj. 
In all cases, we should have Axj = {x^ ~ )• 
We are going to define dj, the coupling strength to probe- j, by 



dj = d{xj)^/Mj 



(71) 
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where d{x) is a real function defined on the decoherence interval. It has dimensions of Energy xLength"'^/^. Similarly, 
the potential of probe-j, Vj, has to be chosen as a continuous function of position of contact, Xj. Let V{x) denote 
this function, i.e., Vj = V{xj). The velocity at probe-j, vj, will then be 

vj = v{xj) = ^2{E-V{xj))/m* . (72) 

Then we will define D{x) function as 

D{x)=^^ , (73) 

and the coeSicients Dj becomes D{xj)y^Axj. For this reason, the function D[x) has the dimensions of Time~^/^. 
Hopefully, we are going to demonstrate with numerical solutions that D{x)^ corresponds to the decoherence rate 
It is natural to define the two functions /l(x) and fR{x) as 

Uix) = D{x)Mx) , fnix) = D{x)ipRix) ■ (74) 

In that case we have /lj = f l{x j) H^x j and fnj = fii{xj)^/Axj. (The functions fL,R{x) have the dimensions 
Length"^/^, but fL,Rj are dimensionless.) 
The r matrix is defined in the usual way as 

Tji = 6je + ^ fRj^ fLj^ = Sje +^fR {xj^ ) /l {xj^ ) y^AxjAxe (75) 
to to 

We are interested in obt ainin g a functional form for the F matrix. Note that in discrete form, F~^ is applied to 

the vectors which have V Ax factors in all of their elements. For this reason, let us investigate the general relation 
Uj = Fji b£ where aj = a{xj),y Axj and bj = h{xj)^ Axj. 

/~Ax~- 

a{xj)^/A^j = b{xj)^/A^j + ^^—^ V fR{x<)f{x>) b{xi)Axe . (76) 

*o e 

Eliminating the common factors in square roots we get a functional equation 

aix) = J V{x;y)b{y)dy , (77) 

where 

T{x;y)=S{x-y) + ^fR{x<)fL{x>) . (78) 
to 

Therefore we are going to define functions tl{x) and tr{x) (defined only on the decoherence interval) by 

fL,Rix) = jT{x;y)TL,Riy)dy ■ (79) 
Using these we have tlj = TL{xj)^/Axj etc. Similarly the inverse of F function can be expressed as 

T-\x;y) = S{x-y)- ^TR{x<)TLix>) , (80) 

where 

td = to-^ TRjfLj = to- TR{x)fL{x)dx . (81) 
i •' 

The refiection amplitudes can also be expressed in the same form. 
The transmission probabilities are 

nj = \TL{xj)fAxj=n{xj)Axj , TRj = \rR{xj)\'Ax,=TR{x,)Axj . (82) 
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It is good that the probabihtics turn out to be proportional to the interval length (Note that the probe-j takes care 
of the decoherence on an interval with length Axj ) . The transmission between two different intervals 

Jd -Id 

is also proportional to both of the lengths of the corresponding intervals. 
Next, note that 

mj = 2Re^^^^^^^^ Axj = M{xj)Axj . (84) 

The matrix elements of P becomes 

Pji = SjiM{xj)Axj - 2 ^fl(^<)^-^('^>) AxjAx£ . (85) 
This matrix looks different from T in the way it contains interval lengths. But still we can define a function form 

P{x-,y) = M{x)d{x-y)-2'^^^^^^^^^^p^ . (86) 
So, if W{xj) denotes the electrostatic potential on probe-j, we have 

nix) = Jp{x;y)^dy . (87) 

P{x; y) also satisfies the equation 

j P{x; y)dy = Tl{x) + Tr{x) . (88) 
Finally, it can be shown that the dimensionless conductance g can be expressed as 

(89) 



= Td + 2 j TL{x)dx-2 J J TL{x)P-\x;y)TL{y)dxdy (90) 

= Td + 2 f f TR{x)p-\x;y)TL{y)dxdy (91) 



where P ^ {x; y) is the inverse of P{x\ y) 

' p-\x;y)P{y;z)dy = 5{x-z) . (92) 



IV. SMALL DECOHERENCE RATE 

In here we will assume that the coupling strength expression d{x) is small, so that we can expand all relevant 
quantities in D{x). Mostly we will be interested in the lowest order term. The functions Jl and fa are of first order 
in g. The T function-matrix is 

T{x;y) = 5{x-y) + ^fR{x<)fL{x>) , r-\x;y) ^ 6{x - y) - ^fn{x<)fL{x>) • (93) 
to to 

From here we get tl « /l and tr « fn where the corrections are of third order. 
The direct transmission amplitude is 



td^to 



J fL{x)fR{x)dx . (94) 
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The direct transmission probability becomes 

(95) 



Note that 



M(x)=2Re^5i^^«2ReM^l^ , (96) 
td to 

which is of second order, as a result we can express as 

Td«|to|'(l- J M{x)dx^ . (97) 

The transmission probability densities to the probes are 

mx)^\fL{x)f , Tn{x) ^ \fR{x)f , (98) 
which are of second order. Therefore, the P matrix-function 

P{x;y) = M{x)S{x-y)-^TR{x<)TLix>) , (99) 

has at least a second order term as the first term and a fourth order term in the last term. For this reason, we might 
need to calculate M{x) to fourth order as well. Let us consider the problem in the following way. Write the matrix 
as P = Pi + P2 where Pi = M{x)S{x — y) and P2 is the remaining term. Inverse of P is 

P-l = Pf 1 - P{^P2Pi^ + P^^PiPi^P^Pi^ (100) 

Since Pj~^ = M{x)~^5{x — y), we have 

T r ^P-V ^T, I . f Tr{x)Tl{x) ^ [ f TRix)P2{x;y)TL{y) ^ ^ , 
TRix)P ix; y)TUy)dxdy = J dx -j j dxdy + • • • , (101) 

where the first term is of second order and the second one is of fourth order. We keep the first term only. For this 
reason, we don't need to calculate the higher order terms in M(x). The result for the dimensionless conductance is 



ItolMl- /M(x)rf.)+2/^«^gMrf, . (102) 



Summary of the steps of a numerical computation 

• A potential V{x) has to be chosen and the solutions (pL,R of the Schrodinger equation at a selected energy E 
have to be obtained. We will use <fiL,R = \/vlV^l,r which are dimensionless. Through the solutions, we also 
obtain the scattering matrix of the "bare" main conductor, the amplitudes tq, t'q and to; but we need only the 
transmission amplitude to- 

• A decoherence interval (from x^ to x^) has to be chosen and a function D{x) has to be defined on this interval. 
D{x) has the dimensions of Length"^/^. It is related to d{x) through the relation D{x) = d{x)/'h^Jv{x)vL- We 
ignore the energy dependence of D and for all energies, E, use the same function. 

• For the calculation, we divide the interval into N subintervals each with length Axj and positioned at 
Xj. We choose N to be large enough so that each subinterval is smaller than the wavelength of solutions (or 
smallest length scales associated with the wavefunctions <fL,R)- 

• We define x 1 column matrices Jlj and Jrj by 

Ilj = D{xj)(fL{xj)y/Axj , fRj = D{xj)ifiR{xj)y/Axj , (103) 
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We construct the T matrix by 



^ji^S.e+^fRjjLj, ■ (104) 



We obtain A'^ x 1 column matrices tlj and Tjij by = F ^f^ and th — T ^ fn 

R 



The direct transmission amphtude is calculated by using td = to — TrIl and the associated probability by 



Td - \tdt 

• The transmission probabilities from the left and right leads to the probes are obtained by T^j = |'''Lj|^ and 
Trj = kfljl^- Also, we find rrij by 

^2Re^5^ . (105) 
td 



We will define a matrix P by 



2 

Pjt = rrijSjg - TfrTRj^TLj^ . (106) 



The dimensionless conductance and the local electrostatic potentials of probes are calculated by 



g - Td + 2T^P-'TL , (107) 



Wl 



^ - [P-'Tl), . (108) 



V. RESULTS AND CONCLUSION 

In this work we have revealed our continuum model for decoherence in ID transport through a mesoscopic wire. 
The dephasing effects in ID transport had been investigated by extending Biittiker dephasing model, which is a 
conceptually simple model to simulate the dephasing effect in ID transport through a mesoscopic system by coupling 
electron reservoir to the conductor. In our model decoherence proceeds at every location such that we coupled 2N 
electron reservoirs to the conductor by 2N channels and we choose N to be large to obtain a continuum case. In the 
reservoirs inelastic events and phase randomization take place. Electrons can go to equilibrium in those channels but 
will eventually return back into the system and at the end, as a result of dephasing, coherence is lost, same as in 
the Biittiker's dephasing model. Our model is more consistent with the prevalent notions of decoherence since the 
placement of the single scatterer in Biittiker's model effects the electron transmission. 

The key point that we have solved in this work is whether extending Biittiker's fictitious probe model can be made 
and give us more reliable data. We apply our model of continuum decoherence for the double barrier case in a one 
dimensional wire at mesoscopic scales and focus on resonant tunneling seen in such devices. 

Incident electrons are described by plane waves. We consider potentials with V{x — oo) = V{x +oo) so that 
= fci?. and vl ~ vji. In this case for (p^ and (pfj can be absorbed into 7, i.e.. 





for 
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—00 


'toe}^'"") 


for 
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+00 


toe-''''') 
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ipR = \fvLVR. = 

= --^= so /jjj- = '^j(p-ii(xj). In this case (ph.R and 7 are dimensionless. 

Electron waves tunnel through the left and right barriers via the quantum well. The potential felt by the electrons 
is depicted in Fig. [2] In the well the electron wave experiences multiple refiections due to the barriers and then 
the wave tunnels out the right barrier. Transfer matrix method is used to calculate the reflection and transmission 
coefficients. The barriers transfer matrices are obtained by matching the wave functions and their derivatives at the 
boundaries. So we had the transmission and reflection amplitudes. Once we get the transmission probability we apply 
our procedure to get the conductance g. 

Fig. [3] shows conductance versus Ep graph for different D values for the double barrier case. As seen in the figure 
the conductance decreases with the increase in decoherence. D=0 case is shown at the top. The peaks seen in the 
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FIG. 2: Double barrier case. 




FIG. 3: Conductance vs Ef graph for different D values. D values are 0, 0.3, 0.5, 0.7, 0.9 from top right to bottom right. 

tunnelling region, where the energies are smaller than Vq = 2.5, are due to resonant transmission. In this region we 
see that decoherence makes the constructive interference of electron waves disappears. After that region we see that 
conductance, i.e. the electron transmission, is suppressed by dephasing. 

Fig. |4] shows conductance versus D graph for Ep ^ Ei = 0.96 and for Ep = E2 ^ 1.41 which is the second 
maximum and second minimum at Conductance vs Ep graph for different D values (Fig. [s]). 

Decoherence mainly prevents wave interference. Depending on whether the interference increase or decrease the 
transmission probability, decoherence may decrease or increase the conductance. So, if constructive interference is 
present in the forward direction decoherence will prevent that and decrease the conductance. Otherwise, if destructive 
interference is effective in the forward direction, then decoherence increases the conductance. But as a rough guide 
we can give the following rule: When the transmission probability is roughly below 0.1, decoherence increases the 
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FIG. 4: Conductance vs D graph for Ep — Ei — 0.96 which is the second maximum at Conductance vs Ep graph for different 
D values(Fig. [3| and for Ep — E2 = 1.41 which is the second minimum in the same Fig. ^ 



conductance. Otherwise, if the transmission probabihty is above 0.1, then decoherence decreases the conductance. 

In summary, we have proposed a model to address the significant dephasing effects in ID transport. And we observe 
that dephasing can dramatically suppress the conductance of a conductor since it effects the transmission probability 
of the electron waves. 
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